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Theory of charmless radiative B decays is reviewed. Existence of uncontrolled 
non-perturbative effects in the inclusive rate at Oicis) is reminded. 



Charmless radiative B decays are the decays B — > X^o charm7, where 
-''^no charm is either a particular hadronic state for an exclusive decay, or just 
any charmless hadronic state in the inclusive case. Such decays are generated 
by tree-level b u transitions with photon radiation, loop-mediated b — » 
d transitions, and loop- mediated b —^ s transitions. Examples of diagrams 
contributing to each of these three types of transitions are presented in Fig. |l|. 
The loop-mediated transitions are known to be very sensitive to new physics, 
e.g. to existence of SUSY particles with masses below 1 TeV. In the Standard 





Figure 1. Examples of Feynman diagrams contributing to B ^ Xno charm7- 



Model, the b 



and b d transitions are CKM-suppressed with respect 



to the b ^ s ones. The relative suppression factors are \Vub/Vtt 



1% and 



|Vtti/Vts| G [2.5%, 5%], respectively. Therefore, to a good approximation, 
the loop-mediated 6 — > s transitions saturate the inclusive K. — > X^o ciiarmSj 
branching ratio. This branching ratio is measured by CLEOEl and ALEPHH 
at the level of around 3 x 10^"*, after the contribution from intermediate ip 
states is subtracted. 
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The latter term in the above equation is the very contribution from interme- 
diate ^p. A lower bound on its numerical size can be found by summing up 
the exclusive branching ratios of charmless radiative decays listed in the 
Particle Data BookEl (which give together around 4%) and multiplying them 
by BR[B Xip] that is close to 1%. It follows that the intermediate ijj con- 
tribution to BR[B Xsj] is not smaller than 4 x 10~^, i.e. it is larger than 
the remainder of this branching ratio. 

The CLEO result is interpreted here as the one with subtracted interme- 
diate ip contribution, even though no subtraction has. been actually made in 
the measurement performed for high-energy photons However, the extrapo- 
lation to lower photon energies has been done with use of a theoretical model 
that did not include intermediate ip. Photons originating from the intermedi- 
ate ip are expected to be softer than most of the other photons in i? ^ Xs^- 
A quantitative estimate of their softness (missing at present) is necessary to 
completely clarify this point. 

The intermediate "ip contribution had to be included in Eq. (^, because 
each Xno charm IS assumcd to be a QCD-eigenstate, while ip is not stable in 
QCD. In other words, the diagrams in Fig. |l], when dressed by an appropriate 
number of gluons, give a contribution to the intermediate ^p channel as well. 
Thus, a separation of the intermediate ip contribution, which may be straight- 
forward on-the experimental side, has to be thought about on the theoretical 



Most of this talk will be devoted to the loop-mediated B Xg^ decay 
that is dominant in i? — > X^o charm7. In order to make a theoretical pre- 
diction, we first need to calculate the perturbative &-quark decay amplitudes 

(v) 

to partonic states Xg and the photon. Later, the perturbative amplitudes 
will enter directly into the expressions for hadronic branching ratios. A lot of 
effort has been devoted in the recent years to calculating the &-quark decay 
amplitudes with better than 10% accuracy. Single gluon corrections (Fig. ||b) 
to the one-loop 6 ^ 57 diagrams (Fig. ^) increase the predicted amplitude 
by around 50%, and the branching ratio by around 100%. This effect is so 

large because the logarithm In — ^ is big and because the one-loop result is 

accidentally quite small - it gives only about ^ of what is naively expected. 
In order to achieve better than 10% accuracy, one needs to include the NLO 
QCD corxections (Figs. ||c and ||d), i.e. non-logarithmici parts of two-loop 
diagramsQ'Ll and logarithmic parts of three- loop diagrams The NLO correc- 
tions further increase the predicted branching ratio by around 20%. Both the 
LO and the NLO calculations include resummation of large logarithms In — W- 
from all orders of the perturbation series. 
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Figure 2. Examples of Feynman diagrams contributing to 6 — > 57 at various orders in the 
renormalization-group-improved perturbation theory. 

Rcsummation of large logarithms as well as further calculation of hadronic 
decay rates is most conveniently performed in the framework of an effective 
theory obtained from the SM by decoupling the heavy electroweak bosons and 
the top quark. The Lagrangian of the effective theory reads 



^ = ^QCDxQED {u, d, s, c, b) + -^VtlVtb Ci{ii)Oi 

^ i=l 



(2) 



where the first term is just the QCDxQED Lagrangian for the light quarks, 
and the second term contains flavour-changing local interactions Oi of either 
4 quarks or 2 quarks and gauge bosons. 



Oi 



{rsT,b)j:,{qm, 



1,2, 

3,4,5,6, 
7, 



\C^{^lb)\ - 1, 

\aifib)\ < 0.07, 

|C8(/x,)| ^0.15. 



(3) 



The symbols Fj and F- in Oi,..., Oq stand for various products of the Dirac 
and colour matrices. The MS-renormalized couplings Cj at the scale /Xb ~ 
nib are known nowadays up to (and including) the following terms in their 
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Once the Wilson coefficients Cidib) are known, the 6-quark decay amplitudes 
are given by Feynman diagrams with single insertions of the flavour-changing 
interactions, i.e. by matrix elements of the operators Oi between the appro- 
priate partonic states. 

For the exclusive decay B — > K*j, one needs to know matrix elements 
of those operators between the relevant hadronic states: {K*j\Oi\B) . There 
have been many attempts to calculate these matrix elements using quark 
models, QCD sum rules, lattice and heavy quark symmetries (see e.g.EHl). 
The history of published predictions is briefly summarized in Fig. ^j. The 
(blue) thin bars and dots are the theoretical predictions, while the (red) thick 
bars denote the CLEO measurements. Many recent theoretical papers on B ^ 
K*^ are not included in the plot because only form-factors are discussed there, 
and no explicit number for the decay rate is given. A general conclusion one 



BR[B K*-f] X 10^ 



I. 



199B 2000 



Figure 3. Brief history of published predictions and experimental results for BR[B K*^]. 



can derive from the B — > K*^ studies is that they can help us in understanding 
non-pcrturbative QCD, but BR[B K*'^] is not a good place to look for 
new physics, given the present experimental and theoretical results for the 
inclusive mode B — > Xgj. However, observation of a sizablCpCP-asymmetry 
in the exclusive mode would be a clear signal of new physics JHI 

For the inclusive decay B Xs^, the theoretical prediction for the 
branching ratio can be made more precise by using the Operator Pmduct 
Expansion (OPE) within the Heavy Quark Effective Theory (HQET).0 We 
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need to calculate 



E 



Cr{X,^\07\B) + C2{X,^\02\B) + ..Y , (5) 



where dots stand for matrix elements of other operators that are numerically 
less important. Let us first look at the "77" term, i.e. the term proportional 
to IC7P in Eq. (||). This term dominates in the perturbative calculation of 
the &-quark decay. In full analogy to the semilcptonic _B-meson decay, we 
relate it via optical theorem to the imaginary part of the elastic forward 
scattering amplitude 



} = ImA (6) 



In this amplitude, we can perform OPE when the photon energies E-y in the 
_B-meson rest frame are far from the endpoint, i.e. when jms — 2E-y\ » 
Aqcd- Most of the photons in B Xg-f have energies close to the endpoint 
^max ^ ijTjj,^ so they do not satisfy this requirement. Thus, at this first step, 
OPE gives us only the tail of the photon spectrum. 




2 E:™^'^ \ Re E^ [GeV] 



Figure 4. Physical cuts in A(Ej) (thick lines on the real axis) and the relevant contour. 

Fortunately, we know analytic properties of the amplitude A when E-y is 
formally treated as complex. We know that the discontinuity of A on the real 
axis is equal to ImA. Thus, if we want to know the integral of ImA from, say, 
1 GeV to the endpoint, we can find it by performing the integral of A around 
the big circle in Fig. ^, where the condition for OPE is always fulfilled. 



'1 GcV -/big circle 



dE^ E" Im A{E^) - 6 dE^ E!" A{E^), n = 0, 1, 2, ... (7) 
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We should better not go with much below 1 GeV. There is no problertLin 
doing so for the "77" term, but there are problems with other operators.l3il3 
Anyway, the region below 1 GeV is hardly accessible experimentally, because 
of the b c background. 

The conclusion at this point is that we can predict the photon spectrum for 
not too small and not too big energies, and moments of the photon spectrum 
from not too small energies to the endpoint. Making such predictions requires 
calculating matrix elements of various local operators among S-meson states. 
Matrix elements of higher-dimensional operators are suppressed by higher 
powers of A/ms- Therefore, we can write a double expansion for the "77" 
term, i.e. we can write 

^ BR[B ^ X,7] 



+ [ Contributions other than the "77" term]. (8) 

Here, the two terms not suppressed by A/ms are simply those already found 
in the perturbative calculation of the 6-quark decay. The term proportional 
to ao2 turns out to give only an around —3% contributionJlj The remaining 
terms in the first two lines of Eq. (||) have stronger suppression factors, which 
makes them negligible. 

However, we need to ask whether a similar expansion can be written for 
the third line of Eq. (||), i.e. for the contributions other than the "77" term. 
The answer to this question is no. These remaining contributions contain, for 
instance, the huge effect from intermediate "0 that has been mentioned in the 
beginning of this talk. The intermediate ip becomes important either due to 
non-perturbative effects or because of big contributions at high orders of per- 
turbation theory that need to be resummed. Most probably, both mechanisms 
are at work. 

The intermediate ip contribution can be just subtracted from both the 
experimental data and the theoretical predictions, using the narrow peak ap- 
proximation as in Eq. (^. Other narrow cc resonances hardly ever decay 
radiatively to charmless states, so similar contributions from them are negli- 
gible. 

Suppose we subtract the intermediate ip contribution. Does the sum 
of the remaining contributions take the form of a power series as in the 
first two lines of Eq. (||), with the perturbatively calculable leading term? 
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It does not.iS However, it is hard to identify any obvious source of a 
big non-perturbative effect in it. Operators containing no charm quark 
are suppressed by their small Wilson coefficients. As far as the operators 
containing the charm quark are concerned, we know that their contribution 
at the leading order in ag can be expressed as a power series 



{B\ V.y ^ \B) = (perturbative 0) + ^ V 6„ , (9) 

O2 O7 mi^^ \mi J 

which can be truncated to the leading|jj— p term, because the coefficients 
bn decrease fast with n. The calculablda'tll n — term makes BR[B ^s7] 
increase by around 3%. 

However, an analysis of non-perturbative effects in the matrix elements 
of Oi and O2 at 0{as) is missing. For instance, 




{B\ ^ . -^ l-B) — ^ono-loop + B^ + Cunknown, (10) 

O2 O2 

where ^ono-ioop stands for the very small (less than 1% in BR) perturbative 
contribution from the gluon bremsstrahlung at one loop, B^ is a part of 
the (huge) intermediate ip contribution, and Cunknown denotes the remaining 
non-perturbative terms. Those remaining terms would not be numerically 
important if they were either suppressed by A/rric^b, or small for other reasons, 
or could be absorbed into the intermediate ip contribution. Unfortunately, I 
am not aware of any sufficiently precise argument that any of these three 
possibilities is realized. 

In the following, I shall assume that one of these three possibilities is 
realized. In such a case, the hadronic decay rate is indeed well-approximated 
by the partonic decay rate, up to small non-perturbative corrections 

T^[i=i V 1 subtracted -0 t^tl v ipcrturbativc NLO 



X [1 + (O(AVmg) ~ 1%) + {0{K^lml) ~ 3%)] . (11) 

The normalization to the semileptonic rate has been used here to cancel un- 
certainties due to , CKM-angles and some of the non-perturbative correc- 
tions. One has to remember that Eq. ( pT| ) becomes a bad approximation for 
i?™* << 1 GeV, and that non-perturbative corrections grow dramaticallyO 
when i;™* > 2 GeV. 
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For Er 



= 1 GeV, Eq. (|T|) gives 
BR[B 



"v 1 subtracted ib 



(3.29 ±0.33) X 10" 



(12) 



where the dominant uncertainties originate from the uncalculated 0{a^g) ef- 
fects and from the ratio mc/mb in the semileptonic decay (around 7% each). 

Unfortunately, E'cut = 1 GeV is not accessible experimentally at present. 
We need some prediction for tiie photon spectrum. The solid line in Fig. ^ 
describes the photon spectrunJla in the region where the theoretical HQET 
prediction is solid. For larger energies, the less solid the line becomes, the less 
solid the prediction is. In the peak region, it is simply an "artist view" of how 
the spectrum could look like. However, its normalization is fixed by Eq. (p^), 
and the_size of the visible K*{%92) peak is adjusted to the value measured by 
CLEO.y 



subtracted ip IQ"^ 



[GeV] 



Figure 5. An "artist view" of ^B_R[B x s'y]'"'^^'^'"'^'"^ 

If we knew the shape function of the B meson exactly, we ccmld make a 
solid prediction for the photon spectrum in the peak region, too.llZI Unfortu- 
nately, only models for the shape function are available at present. Therefore, 
an optimal scenario for a comparison between theory and experiment would 
be measuring the photon spectrum above ^ 2 GeV without relying on any the- 
oretical prediction for its shape, and then extrapolating in a simple manner to 
the predicted snectrum below 2 GeV. The measured spectrum above 2 GeV 
would providcLj important information for extracting Vub from B — s- X^lv- 
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The decay B Xd 7 is theoretically more difficult than i? — > 7, 
because diagrams with up-quark loops are no longer CKM-suppressed with 
respect to the remaining ones. Thus, the theoretical accuracy is at best ±30%, 
even for fixed Wolfenstein parameters p and 77. 

When the noiiiperturbative effects in up-quark loops are assumed to be 
small, one obtainsllj 

\{BR[B -> 7] + BR[B ^ Xd 7]} ^ 2.43[(1 - pf + f - 0.35(1 - p) 

+0.07] X 10-^ = 1.61 X 10-^ (for p = 0.11 and f] = 0.32), (13) 

while the direct CP-asymmetries range from 7% to 35%. Estimates for the ex- 
clusive channels are: BR[B /9±7] e [1,4] x 10"^ and BR[B (p",w)7] G 
[0.5,2] X 10-6. 

The present experimental results for B —> X^^ already place severe con- 
straints on extensions of the SM, like 2HDM, MSSM, LR- models etc. Theoret- 
ical predictions for exotic contributions have been recently calculated at NLO 
in many extensions of the SM.Ej These NLO effects are important only when 
the exotic effects are large, but tend to cancel among each other and/or the 
SM contribution, so that the present experimental constraints are satisfied. 

The CP-asymmetp£ in _B — > K*^ is very small in the SM, but could be 
significantly enhancedliil in such extensions of the SM, in which the flavour- 
changing interactions of the right-handed s-quark are not suppressed tw 
ms|Mw^ e.g. in the left-right symmetric models. The present CLEO boundll3 
on the CP-asymmetry places important constraints on such models. Interest- 
ing information on physics beyond the SM canJpe obtained from the CP- 
asymmetry in the inclusive B — > Xs^ mode, tooHd 

To conclude: 

• The present theoretical prediction for BR[B Xg'^] in the SM agrees 
very well with the measurements of CLEO and ALEPH. However, an 
analysis of non-perturbative effects at order 0{as) is necessary in order 
to make sure that the theoretical uncertainties are indeed around 10%. 

• Future measurements of BR[B Xgj] should rely as little as possible 
on theoretical predictions for the precise shape of the photon spectrum 
above - 2 GeV. 
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